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ABSTRACT 

We combine newly measured rotation velocities, velocity dispersions, and stellar masses to con- 
struct stellar mass Tully-Fisher relations (M^TFRs) for 544 galaxies with strong emission lines at 
0.1 < z < 1.2 from the All Wavelength Extended Groth Strip International Survey (AEGIS) and 
the Deep Extragalactic Evolutionary Probe 2 Survey (DEEP2). The conventional M*TFR using only 
rotation velocity (K-ot) shows large scatter (~ 1.5 dex in velocity). The scatter and residuals are 
correlated with morphology in the sense that disturbed, compact, and major merger galaxies have 
lower velocities for their masses. We construct an M*TFR using the kinematic estimator 5*0.5 which is 
defined as ^0.5V^^f- + cr^ and accounts for disordered or non-circular motions through the gas velocity 
dispersion (cTg). The new M*TFR, termed Sq.s/M^TFR, is remarkably tight over 0.1 < z < 1.2 with 
no detectable evolution of its intercept or slope with redshift. The average best fit relation has 0.47 
dex scatter in stellar mass, corresponding to '--^ 1.2 'magnitudes,' assuming a constant mass-to-light 
ratio. Interestingly, the ^'o.s/M^.TFR is consistent with the absorption-line based stellar mass Faber- 
Jackson relation for nearby elliptical galaxies in terms of slope and intercept, which might suggest a 
physical connection between the two relations. 

Subject headings: galaxies: evolution - galaxies: general - galaxies: high-redshift -galaxies: interac- 
tions - galaxies: kinematics and dynamics - galaxies: spiral 



1. INTRODUCTION 

The Tully-Fisher relation (TFR) between the luminos- 
ity of galaxies and their rotation velocity is a fundamen- 
tal scaling relation that constrains galaxy formation and 
evolu tion models (e.g., iDalcanton. Spergel. fc Summers! 
Il997l and references therein). Recent local stu dies fo- 
cus on stellar or bary o nic masses (e.g., [Be ll &: de Jonel 
120011: iMcGaughl [200a iPizagno et all l2005f) . which are 
easier to model than luminosity. Distant stellar-mass 
TFRs (M»TFRs) have also been studied (Conselice et 
al. 2005; Flores et al. 2006). Neither found evidence for 
evolution to redshifts z ~ 1 or 0.6, respectively, but each 
pruned their data to exclude morphologically odd galax- 
ies. This letter presents a study of the M^TFR to redshift 
z = 1.2 using a sampl e of 544 galax i es with strong emis- 
sion lines. Following IWeiner et al.l ()2006bf l. we exclude 
few galaxies and adopt a new velocity measure that in- 
cludes both velocity gradients (usually rotation) and ve- 
locity spread (dispersion) . We adopt a ACDM cosmology 
with h = 0.7, = 0.3, and flA = 0.7. 



Galaxy kinematics and morphologies at z ~ 1 are pos- 
sible with high resolution spectra and images. To track 
different morphologies and redshifts, samples should be 
large enough to be subdivided. The All- Wavelength Ex- 
tended Groth Strip International Survey (AEGIS) pro- 
vides such data. Single-orbit images were taken by the 
Hubble Space Telescope Advanced Camera for Surveys 
(HST/ACS) in the F606W {V) and F814W (/) band- 
passes, and spectra were obtained by the Deep Ex- 
tragalactic Evolutiona r y Probe 2 S urvey (DEEP2); see 
details in I Davis et"aI1 (|2003l l2006f ) . Obs ervations and 
reducti on of kinematics are described in IWeiner et al.l 
(|2006a^ for data similar to DEEP2. Kinematics are mea- 
sured from the strongest emission line(s) among Ha, H/3, 
[OH] A3727, and [O HI] A5007. The spectral resolution is 
FWHM==1.4 A, or 56 km at z = 1 (Gaussian cr = 24 
km s^^). Emission line widths can be measured down to 
~ 0.6 of the cr or ~ 15 km s~^, and velocity centroids 
for rotation curves are good to 0.4 of tr or ~ 10 km s^^. 
Rotation velocities can be detected down to ^ 5 km s~^, 
if present. 
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3. SAMPLE SELECTION 

We chose a sample of galaxies that span a wide range in 
redshift, stellar mass, and morphology. The DEEP2 sur- 
vey had magnitude limits of Rab ~ 24.1, bright enough 
to yield good HST/ACS imaging. A redshift limit of 
z = 1.2 e nsures that galaxies of all colors on the "blue se- 
quence (iWillmer et al.l[2006h are well-sampled. To yield 
good quality kinematics, emission lines were required to 
have integrated intensities > 1500 e~/A in the summed 
one-dimensional spectrum, and spectrographic slits had 
to be aligned to within 40° of the major a xes of the 
galaxies (see Fig. 13 of IWeiner et al.l [2006a[ ). Also, to 
reduce the effects of dust, edge-on systems with i > 70 
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were excluded, and to allow a reliable rotation measure, 
nearly face on galaxies with i < 30 were removed. A total 
of 14 galaxies were excluded because they had emission 
lines with disturbed morphologies or that had contam- 
ination by emission from a nearby galax y on the sky. 
After addition al cosmetic cuts were made (|Kassin et al.l 
lin preparationf ) ■ the final sample became 544 galaxies. 
Our final sample is primarily selected on emission line 
strength. 

4. DATA REDUCTION AND ANALYSIS 

4.1. Photometric Parameters 

In order to interpret dynamical measurements made 
from the spectra, accurate position angles and elliptici- 
ties of galaxies must be obtained from HST imaging. To 
measure these para meters, the S Extractor galaxy pho- 
tometry soft ware (iBertin fc Arno uts 1996) is used, as 
discussed in iLotz et al.l ( 2006f ). To obtain quantitative 
m orphologies for galaxies, we use t he Gini/M2o method 
of lLotz. Primack. fc Madaul ()2004D . This system of ob- 
jective structural parameters has been shown to reli- 
ably divide galaxies into bins analogous to Hubble types; 
namely, early type spirals, late type spirals and irregu- 
lars, and major mergers. Galaxies are classified in the V 
image for z < 0.6 and in / for z > 0.6. Such quantitative 
morphologies are obtainable only for galaxies with Pet- 
rosian radii > 0.3", images that have a signal-to-noise 
ratio of > 2.5, and that do not fall near the edge of 
the HST/ ACS CCD chip. Galaxies without quantitative 
morphologies are not removed from the sample for com- 
pleteness sake, and do not differ statistically from the 
rest of the sample. All photometric parameters, includ- 
ing quanti tative morphologie s, will appear in a forthcom- 
ing paper, iLotz et al] (|2007( ). In a separate exercise, we 
visually classified all galaxies and found examples that 
were fairly normal according to Gini/M2o, but looked 
to our eye to be more disturbed or compact; we discuss 
these galaxies separately. 

4.2. Spectroscopic Measurements 

Spatially extended emission lines in the galaxy spectra 
are used to measure gas rotation and dispersion profiles. 
Since the spatial extent of the line emission is at most 
only a few times the seeing (typically 0.7" FWHM), the 
effect of seeing must be modeled. We do this with th e 
ROTCURVE fitting procedure of IWeiner et"aI1 (|2006aD . 
The kinematic model we fit has two parameters: the ve- 
locity on the flat part of the rotation velocity profile, Vrot , 
and the velocity dispersion from a spatially resolved fit to 
the two-dimensional spectra, ag . The best- fit Vrot values 
are corrected for inclination with the measured galaxy 
ellipticities. For galaxies without a disk-like geometry 
(based on visual inspection of HST/ ACS images), Vrot is 
not inclination-corrected. These galaxies are flagged in 
the following analysis, but fits do not change if all galax- 
ies are inclination-corrected. The ag we measure does 
not have to correspond to a literal gas velocity disper- 
sion like that of stars in an elliptical galaxy. Instead, for 
(Tg < 20 km s~^, <7g likely measures the relative motions 
of individual H II regions, as in spiral arms or a thick 
disk; for ag > 20 km s~^, ag likely represents an effec- 
tive velocity dispersion caused by the blurring of velocity 
gradients on scales at or below the seeing limit that may 



not even have a preferred plane. An uncertainty of 0.1 
dex is adopted in both ag and the inclination-corrected 
Vrot to account for random errors and the dependence of 
the model parameters on the assumed seeing and scale 
radius of the rotation curve. Results do not change sig- 
nificantly if uncertainties of 0.2 are chosen, as discussed 
in §6. 

5. STELLAR MASS ESTIMATES 

A few methods of estimating the stellar masses (M*) 
of galaxies have been invest igated in a separate study 
(|Kassin et al.l lin preparation!). These m ethods are: a 
color-M/L relation from lBell et al.] (I2003D for ]/ and 
M/Lb, a color-M/L relation from iBell et~al] (|2003D for 
B — V and M/Lh, M/Lh = 1, broad-band spectral en- 
ergy d istribution (SED) fits of lBundy. Ellis, fc Conselici 
(|2005D which incorporate observed o ptical and K-h&Yid 
data, and the method of iLin et al.l (I2006D . We choose 
to adopt the method of lLin et al.l ( 20061 ). These authors 
calibrated rest-frame UBV photometry and redshifts to 
the fuU SED fits of Bundy ct al. (2005) to obtain es- 
timates for very blue galaxies that are not detected at K , 
which Bundy et al. (2005) require to derive a M*. While 
the evolution of SEDs with redshift is included and er- 
rors of 0.2 dex are assumed, the results of this paper, as 
discussed next, change insignificantly if we do not take 
into account evolution or inc rease the errors to 0.3 dex 
(|Kannappan fc Gawiseij[2007[ ). 

6. 5o 5 AS A TRACER OF GALAXY-DARK HALO 
POTENTIAL WELLS 

The quantity S*^ = KVrot + o'gi where K is a, con- 
stant, combines dynamical support from ordered mo- 
tion w ith that from disordered motions (|Weiner et al.l 
l2006al) . For a spherically symmetric tracer population 
with isotropic velocity dispersion a.nd de nsity oc r~", 
ag = Vrotl\foL (jBinnev fc Tremaind ll987L §4.4), where 
K = 1/a. Therefore, if galaxies are virialized systems, 
and these assumptions about the tracer population are at 
least approximately correct, then Sk should at least ap- 
proximately trace the global a of galaxy-halo systems. In 
this case, the tracer population is the gas producing the 
emission lines, which we assume follows the M*. For disk 
galaxies with exponential M* profiles, a 2 — 3 brack- 
ets a range of reasonable values (§5.2 of IWeiner et al.l 
l2006al) . A satisfactory approximation is adopted for all 
morphologies, K = if K = ^ is adopted, the overall 
results do not change. We henceforth use S'0.5 for our 
analysis. 

7. RESULTS 

The top panel of Figure 1 shows Vrot versus (the 
Vrot/M^,TFR) for a range of redshifts. The majority of 
spirals (green triangles and blue squares not outlined 
in black) form a clea r ridge-line that comp ares well to 
a local relation from iBell fc de Jond ()2001l ) , consisting 
of well-ordered spirals , and to a z ~ 1 relation from 
IConselice et al.l ()2005l ). who included only "elongated 
disks." In contrast, almost all of the other classes (red 
circles and symbols outlined in black) lie to lower Vrot, 
causing a large scatter ('^ 1.5 dex in velocity). 

This morphological dependence of TuUy-Fisher scat- 
ter is well known. Among local studies, e.g., larger 
scatter is found for close pairs with kinematic disorders 
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Fig. 1. — The 6*0.5 and Vrot stellar-mass Tully-Fisher relations (M^TFRs) for 0.1 < z < 1.2 in bins of z. Galaxies 
are plotted as different symbols according to their quantitative morphologies, and if a morphology is unobtainable, it 
is placed into the 'none' category. Symbols are outlined in black if the galaxies were considered to be disturbed or 
compact via a separate visual classification. If Vrot is not inclination corrected, open symbols are used; all galaxies in 
the 'none' category are inclination-corrected as their symbol does not have an open form. Using 5*0.5, which combines 
the dynamical support from ordered motion with that from disordered motions, results in an M^TFR that is much 
tighter and non-evolving to z — 1.2. For 5*0.5, a fit to data in the lowest z bin is plotted as a solid line that is repeated 
in the other z bins as dashed lines. Fits are labeled in each box, where a is the intercept at IO^^Mq, b the slope, c the 
intrinsic scatter in 5o.5, and r ms t he rms scatter in 5o.5. In the lowest z bin for Vrot, the M*TFRs from local and z ~ 1 
studies (jBell fc de JonEj[2001l and lConsehce et ani2005l respectively) are plotted as two dot-dashed an d dotted lines, 
respe ctively, to delineate rms scatter. In the lowest z bin for S*o.5, the local Faber- Jackson relation (|Gallazzi et al.l 
I2OO6I ) is plotted as two dotted lines to delineate rms scatter. All relations from the literature have been converted to a 
Chabrier initial mass function when necessary. An estimate of the S*o.5/-Mf,arj/onTFR (§9) is plotted as a long dashed 
blue line in the lowest z bin. 
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(iBarton et al.ll2001[) and when p eculiar galaxies are in- 
cluded (jKannappan et al.1 l2002f ). At higher redshifts, 
larger scatter is found for z ^ 0.5 galaxies with "com- 
plex kinematics" (jFlores et al.| [2006]. Moreover, ~ 25% 
of galaxies found with kinematics unrelated to rotation 
are e xcluded from TuUy - Fishe r studies at 0.25 < z < 
0.45 (ISimard fc Pritchetl 119981) and at 0.1 < z < 1.0 
(jBohm et al.l 120041) . In summary, becaus e of morpho- 
logical pruning, no study except that of iWeiner et al.l 
(|2006bf ) incorporates the population of TuUy-Fisher out- 
liers that we detect in a measurement of Tully-Fisher 
evolution, and the conclusions of such studies will only 
apply for the subsample of galaxies that are morpholog- 
ically well-ordered. 

For M* > 1O^°M0, the percentage of galaxies that 
scatter to low Vrot from the Tully-Fisher ridge-line (de- 
fined as the left-most dotted line in Figure 1) increases 
with redshift (18%, 35%, 42% and 62% for t he four red- 
shift bins plotted). Similarly. IWeiner et al.l ()2006al ) find 
bright galaxies with Cg > Yrot at z > 0.5 that are rare 
at low redshift. 

The bottom panel of Fig. 1 uses 5*0.5 instead of Vrot in 
the M:,TFR. The resulting relation has scatter that is not 
much greater than our measurement uncertainties and 
does not evolve significantly with z.^ Evidently, galaxies 
with low Vrot compared to the K-oi/M,TFR ridge- line 
have a significant ag component which causes them to 
lie on the Sq.s/M^TFR. Lines are fit to the S'o.s/M^TFR 
in the four redshift bins in Fig. 1 with a maximum like- 
lihood method detailed in ' Weiner et all (|2006b( ). We fit 
logio 5*0.5 — a + felogio M, and all excess scatter is allo- 
cated to the 5o.5-coordinate. The fits do not vary signif- 
icantly for the different redshift bins. Average values for 
the slope, intercept at 10^° Mq, and intrinsic scatter over 
all redshifts are 0.34 ± 0.05, 1.89 ± 0.03, and 0.10 dex in 
iSo.s- The average rms scatter is 0.16 dex in 5o.5, which 
corresponds to 0.47 dex in M^, (~ 1.2 magnitudes). The 
fitted relations are not very sensitive to the error esti- 
mates for 5o.5 and A/,; changing either error estimate by 
±0.1 dex changes fit intercepts within the 1-a error, and 
slopes by not more than ±0.1 dex. However, our limited 
range at higher redshift makes slope measurements 
more uncertain than the formal fit errors. Furthermore, 
different estimates of M, can yield a range of ~ 0.2 dex 
in intercept, and estimators with non-evolving SEDs 
can cause the 5o.5/M,TFR to show ^ 0.1 dex evolu- 
tion in intercept with redshift, but this reflects problems 
of these estimators rather than actual evolution in the 
5o.5/M,TFR. 

Since no Tully-Fisher sample at low or high redshift has 
as much scatter as our K-ot/M,TFR due to morphologi- 
cal pruning, it is rather remarkable that the 5o.5/M,TFR 
scatter for the same sample of galaxies manages to ap- 
proach as close as it does to the scatters found for the 
Frot/M^TFRs for galaxies with re latively undisturbe d 
morphologies (~ 0.6 dex in M* for Conselice et al]|2005[ ). 

* One would expect a small (~ 0.05 dex) decrease in intercept 
from z = 1 to 2 = due simply to the predicted increase in M» 
of ~ 50% for a typical spiral galaxy {z formation = 5 and ^.n ex- 
ponentially decreasing star formation rate with t = 7 Gyr), but 
it would not be detectable given our errors. We do not attempt 
to determine K by minimizing the scatter in the Sj^/AfjTFR be- 
cause K is sensitive to the morphological mix of the sample and 
the distribution of measurement errors. 



In addition. Fig. 1 shows that the Faber- Jackson re- 
lation for lo w redshift (0.0 5 < z < 0.22) early type 
galaxies from^ Gallazzi et aTl (|2006[ ) is consistent with the 
'S'o.s/M^TFR in terms of slope and intercept. 

8. DISCUSSION 

If those galaxies with Af, > 10^" Mq that are off 
the Vrot/M^,TFR ridge-line with low Vrot at higher red- 
shift eventually form into well-ordered disk galaxies on 
the ridge-line, this trend implies that more massive disk 
galaxies may be moving onto the K-ot/M^TFR with time 
as their stars and gas "settle" into more circular /planar- 
dominated orbits. The possible mechanisms behind this 
settling can be circularization of gas orbits, gaseous dissi- 
pation to the disk plane, and/or mergers ending with the 
growth of a gaseous disk. The scatter for Af, < IO^^Mq, 
which we can observe at lower redshift, is also large, sug- 
gesting that lower Af* galaxies may currently be settling. 

A scenario of galaxy formation that is consistent with 
our results is one that begins with matter assembling in 
dark halo potential wells with random orbital kinemat- 
ics. The baryonic components form proto-disks that are 
initially supported by a combination of ordered and ran- 
dom motions. The material in these proto-disks has been 
settling since then, unless they undergo major mergers. 
Over the last half or so of their lives, proto-disks and 
their descendants lie on nearly the same S'o.s/M^TFR. 
At first glance, the similarity between the Faber- Jackson 
relation and the S'q.s/Mh.TFR in slope and intercept is 
reassuring, given that merging proto-disks on parabolic 
binary orbits should put the merger pro ducts close to the 
Faber- Jackson relation in Fig. 1 (e.g., iRobertson et al.l 
I2006D . This suggests that the origin of the Faber- Jackson 
relation could in fact be the pre-existing 5o.5/M,TFR 
for disk-like galaxies, if merging is the process that cre- 
ates ellipticals. However, additional studies are needed 
to confirm this speculative scenario. 

9. RELATION TO HALO KINEMATICS 

If the framework we adopt is correct, halo kine- 
matic properties should be related to measurable galaxy 
kinematics. In particular, N-body models show that 
the slope of the relation between halo mass (Mg) and 
halo a (as), both measured within the scale radius 
(iNavarro. Frenk. fc Whit^ll997[ ). is 0.33: login o-^ = 0.33 
logio Ms + constant. Remarkably, this slope is nearly 
consistent with both galaxy relations. In addition, nei- 
ther the halo as — Ms relation, nor Ms or as individually, 
is pre dicted to evolve to w ithin ^ 10% from z = 1.2 to 
now (jPiemand et al.l l2007t see IWechsler etldl |2002| for 
evolution in terms of virial quantities). So, if 5*0.5 is lin- 
early related to halo ag and M^:/Ms =constant, then the 
5o.5/M*TFR should also evolve little from z = 1.2, in 
agreement with observations. Furthermore, that neither 
Ms nor as individually evolves might imply that absolute 
changes in inner baryon structures are also small. 

However, when comparing galaxies to halos, there are 
at least three major complications. First, since Af* is 
subject to the star-formation history of a galaxy, it is not 
as fundamental a quantity as baryonic mass (Mharyon)- 
Locally, low-Af* galaxies have a greate r gas-to-star frac - 
tion than more massive galaxies (e.g.. lMcGaughl l2005fl. 
so the low-redshift So.s/Mhar^onTFR ought to have a 
steeper slope than the 5o.5/M*TFR. Indeed, if we as- 
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sume that Mgas = 0.56Af^ +3.96 (roughly cons istent with 
data from lMcGaughll2005l using relations from lBell et al. 
2003' to obtain M, , which we later convert to a lChabrier 
2003 initial mass function), then the S'o.s/M^TFR slope 
in the lowest redshift bin steepens to ~ 0.39. Secondly, 
galaxies should convert gas into stars over time, causing 
evolution in M^/M^. The third issue is the complicated 
relation between S'0.5 and ct^, which is expected to be 
affected by baryon dissipation and the resp onse of the 
dark halo to it fe.g.. iBlumenthal et al.l lT986'), neither of 
which are completely understood. With a better under- 
standing or measurement of these factors and their time 
evolution, the relationship between galaxy and halo seal- 
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